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In this article I review recent progress towards the analytical calculation of massless 6-point amplitudes. A way 
to organize such calculations is sketched and results for scattering amplitudes in the Yukawa model are presented. 


1. Introduction 

The high energy experiments of the next ten 
years are hadron colliders operating at the multi- 
TeV scale. Because of the high center-of-mass 
energies at the Tevatron and the LHC, the final 
states will look very complex. One may say that 
we are entering the decade of mzdfi-particle/jet 
physics. 

The expected jet rates at the LHC are huge. A 
leading order analysis is not adequate to make a 
detailed prediction for these cross sections. To ad¬ 
dress the problems with a leading order descrip¬ 
tion three points can be made: 

1) Large scale dependence: An N-jet cross sec¬ 
tion behaves like a s (/i) N , means scale uncertain¬ 
ties are amplified with the number of jets. 

2) Peripheral phase space regions: Partonic 
cross sections are very sensitive to next-to- 
leading order effects, when severe cuts are ap¬ 
plied, which is typical for background processes. 

3) Jet structure: The more information about 
the matrix elements is known the better the in¬ 
ternal structure of a jet can be described. 

It is important to note that the ability to de¬ 
tect New Physics at the LHC and Tevatron cru¬ 
cially depends on the understanding of the corre¬ 
sponding Standard Model backgrounds. All that 
is sufficient motivation to treat the generic multi¬ 
particle/jet final states at the next-to-leading or¬ 
der (NLO) level. 

Presently our computational skills for describ¬ 
ing N-jet production at NLO end already at N=3. 
The pioneering work to calculate the amplitudes 
for two-jet processes at NLO in hadronic colli¬ 
sions was accomplished by Ellis and Sexton in 


1986 [ 0. The partonic amplitudes relevant for 3 
jet production were provided by Bern, Dixon and 
Kosower [j|] and Kunszt, Signer, Trocsanyi [| in 
1993/94. New technology was invented for these 
calculations borrowed partly from string theory. 
Further, supersymmetry relations were exploited 
and helicity methods were applied at the NLO 
level [ ||. Other needed ingredients were reduc¬ 
tion formulas for Feynman parameter integrals 
with nontrivial numerators and representations of 
scalar 5-point functions (pentagon integrals) [ . 

For QCD the step to NLO 6-point amplitudes 
has not been taken yet, albeit its phenomenolog¬ 
ical relevance for collider physics. 


2. Scalar scattering in the Yukawa model 


To understand the organization of hexagon am¬ 
plitudes and their efficient computation it is nec¬ 
essary to investigate examples which do not ex¬ 
hibit the full complexity of a generic gauge the¬ 
ory amplitude. An amplitude which contains only 
scalars as external particles is likely to be a good 
candidate. Let us consider 6-scalar scattering in 
the Yukawa model. 
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We restrict ourselves to the massless case and 
the fermion loop contribution. The correspond¬ 
ing amplitude is given by 

r yukbl^2,P3,P4,P5,P6] = 

_gG 
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Each permutation corresponds to a single Feyn¬ 
man diagram. The amplitude for the trivial per- 
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mutation is given by 

■A(pi,P2,P3,P4,P5,P6) = 

d n k tr(/ji fa fa fa fa fa) /o\ 

iW 2 qfqhhhhe 

where qj = k — Tj = k — p\ — ... — pj. Work¬ 
ing out the spinor trace one ends up with scalar 
products of loop and external momenta. One can 
show that all scalar products which contain loop 
momenta can be expressed as inverse propagators 
which can be canceled directly. In this way one 
is getting rid of all loop momenta in the numera¬ 
tor, or in other words: The Yukawa amplitude for 
multi-scalar scattering has a trivial tensor struc¬ 
ture. A given graph of the scalar hexagon ampli¬ 
tude decays immediately into scalar 3,4,5,6-point 
integrals. 

The representation of the amplitude in terms 
of these scalar integrals is not useful as spurious 
infrared poles are present which have to cancel 
in the full amplitude. One can show that none 
of the graphs is actually infrared divergent. The 
tool to make the infrared cancellations manifest 
are scalar reduction formulas. The N-point n- 
dimensional scalar integral integral in Feynman 
parameter space is given by 

I^( Pl ,...,p N ) = (-l) N T(N-n/2) (4) 

jN ^(1 ~ J2l=l x l) 

(x-S-x/ 2) w -"/2 
and the following reduction formula holds 
N 

In = 

k =1 

N 

— (TV — n - 1) E Bk In +2 (5) 

fc =i 

S is the matrix which carries the kinematic in¬ 
formation. Its entries are Mandelstam variables 
Sij... = (Pi + Pj + ■ ■ -) 2 formed by the external 
momenta pi, ■ ■ ■ ,Pn- In-i k stands for the fcth 
pinch integral of I ft, means one leaves out the 
fcth propagator in the momentum representation. 
The reduction coefficients Bk are the solutions of 
a simple linear equation XEi SkiBi = — 1. They 
will play a special role in what follows. They 



can be written very compactly in terms of spinor 
traces. For example 

tr(123456)tr(3456) — 2s34S45S56tr(6123) 

' 4 S 12 S 23 S 34 S 45 S 56 S 61 - tr(123456) 2 

and the other B’s are obtained by cyclic per¬ 
mutations. Here tr(12...) is a shorthand for 
tr( A fa ■ ■ ■), for details see [ 0. 

Applying the reduction formulas allows to ex¬ 
press each Feynman diagram in terms of triangle 
graphs with 0,1 or 2 light-like external momenta 
and 6-dimensional box functions. Only the tri¬ 
angle graphs with light-like legs carry infrared 
poles. As each individual Feynman diagram is in¬ 
frared finite, these triangle graphs have to cancel. 
The cancellations of the infrared divergent trian¬ 
gles gave rise to very interesting insights. There 
are a number of nonlinear constraints between 
Mandelstam variables which are represented lin¬ 
early in terms of the reduction coefficients, e.g. 
one finds 


tr(3456) + tr(123456)Hi + 2 S 34 S 45 S 56 -B 4 = 0 ( 6 ) 

This made the whole structure of the calculation 
very transparent and an astonishingly compact 
result for the Yukawa 6 -scalar amplitude could 
be derived [ fjj . It was found 

r yuk[^1^2,^3,P4,P5,P6] = 

^ ^ G(PlTl 5 P-K2 > Pn3 J Pi** 5 P^5 5 P^Q ) (^) 
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with 


G{Pl,P2,P3,P4,P5,Pe) = g I 3 (Pl2,P34,P56) 

[tr(6123) — 2 s 6 i(si 23 — S12)] 
S123S345 — S12S45 
[tr(1234) — 2534(5123 — S23)] ^ 

S 234 S 123 ~ S 23 S 56 
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+ log(^) log(^)' 

\ s 234/ V 'S.56 J J 

The result consists only of triangle functions and 
log log terms. Apart from dilogarithms hidden 
in the finite triangle function no dilogarithms are 
present in the result. This can be understood by 
the fact that the dilogarithms which are present 
at intermediate steps can be related to single pole 
contributions from triangle graphs with one light- 
like leg and is infrared divergent. As the infrared 
divergence vanishes, the dilogarithnrs also have to 
vanish. The remarkably simple form of the ampli¬ 
tude could be achieved because it was understood 
that the reduction coefficients Bj are the natural 
building blocks of the amplitude. This lead to 
the conjecture that massless hexagon amplitudes 
in general gauge theories can be represented most 
compactly as some linear combination of reduc¬ 
tion coefficients. Still the conjecture has to be if 
not proven then at least confirmed by other calcu¬ 
lations which exhibit a non-trivial tensor struc¬ 
ture. 

3 . 77 —> 4 scalars in the Yukawa model 

To study an amplitude with a non-trivial ten¬ 
sor structure we considered the gauged version of 
Eq. ([!]). For example the fermion loop contribu¬ 
tion to the 

7(Pib(P2) -» <t>{-P3)^>{-Pi)<t>{-Ph)^{-P6) (8) 

amplitude is given by a sum over permutations of 
Feynman diagrams which differ only by the po¬ 
sition of the photons and scalars at the fermion 
loop. 

= g 4e2 juA,a 3 = __ 9 ^_ 

(47r) n / 2 (47r) n / 2 

X 0 ^ ^ 1 2 (Pni 1 PtT2 5 Pn3 5 P'KA 5 P^5 ’ P^6 ) (^) 

Tr&Se 

with, e.g. 

A Xl (pi, Pn 2 , Pn 3 , p T4 , Ptt 5 ,Pn 6 ) = 
d n k tr(/ Al /n /t\ 2 /12 fa fa fa fa) 

:W 2 qhhhhlQe 

Here A j are the helicities of the photons, qj = 
k — Tj = k — pi — ... — pj are the propagator 


momenta. The amplitude exhibits an S 4 Bose 
symmetry, as the scalars are freely interchange¬ 
able. For a given ordering of the scalars one finds 
20 different Feynman diagrams corresponding to 
the different possibilities to attach the two pho¬ 
tons. This set forms a gauge invariant sub-set of 
all Feynman diagrams. It is sufficient to consider 
this sub-set and performing the summation over 
remaining permutations explicitly. 

Gauge and Lorentz structure 

To cope with the combinatorial complexity of 
hexagon amplitudes it is useful to split the prob¬ 
lem into irreducible pieces as far as possible. The 
guidelines are the analyticity and Lorentz struc¬ 
ture of the amplitude together with gauge invari¬ 
ance. The external momenta pi,... ,pe obey mo¬ 
mentum conservation. To make this explicit we 
introduce the five vectors rj = p\ ... + pj , 
j = 1,...,5. We assume for the moment that 
the kinematics is non-exceptional, i.e. the 5 vec¬ 
tors define a 5-dimensional subspace of the n- 
dimensional space. The scattering tensor M has 
the following general decomposition 

M ^ 2 = Ag^+ B hhr£r% (10) 

31 J2=l 

This representation contains 26 coefficients. The 
transversality conditions £1 • pi = 0, £2 ■ P 2 = 0 
reduce the number of relevant coefficients to 17. 
Further, by gauge invariance 

M PlE2 = M £lP2 = 0 (11) 

which leads to 7 independent linear relations be¬ 
tween the tensor coefficients. They may be used 
to express all 17 coefficients in terms of 10 un¬ 
constrained coefficients, e.g. A, Bij, i,j = 3,4,5. 
It follows that one has to calculate only these 10 
independent coefficients. Note that gauge invari¬ 
ance is completely exploited in this way. On the 
other hand, calculating more coefficients and con¬ 
firming the Ward identities explicitly provides a 
test of the calculation. By using helicity methods 
directly one also sees that 10 tensor coefficients 
are sufficient to determine the full amplitude. 

The problem is now reduced to the computa¬ 
tion of 10 tensor coefficients for 20 6-point Feyn¬ 
man diagrams. 
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Projection techniques 

The tensor coefficients in Eq. (10) can be de¬ 
termined by acting with projection operators onto 
each Feynman diagram. This may be achieved by 
defining 


f dnk k _ r (j) 


( 21 ) 


where the parameter integrals with non-trivial 
numerators are defined by [ || 
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( 22 ) 


7+ = 

fc=i 


(13) 


if is the inverse of the 5x5 Gram matrix defined 
by Gjfc = 2rj ■ Tk- Note that the inverse of G 
is only defined if n ^ 4. This will not lead to 
a problem as one can cancel the inverse matrices 
before the limit n —> 4 has to be taken. 

One may easily check the following relations for 
the tensors and 7Z P : 


V pp T > p 1 ' 

= 


(14) 

VWj 

= 0 


(15) 

pi* 

= trip) = tr{V ■ 

■V) = 1 

(16) 

v^n) 

= 0 


(17) 


= fijk 


(18) 
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II 


(19) 

V is a projector onto the (n- 

- 5) dimensional 

sub- 


space perpendicular to the 5 dimensional space 


spanned by the vectors r ? -. TZj^ is the dual vector 
to rj relative to this (n — 5)-dimensional space. 
Here and in the following we use summation con¬ 
ventions for Lorentz and vector indices. 

To determine the tensor coefficients, A and Bki, 
we define the following linear operators 


A(M) = 2 M^ 2 




- 2 H k i Vfi 


AR 


The action of dual vectors 1Z P and projectors 
V^ v on Lorentz indices induces a natural mapping 
from the momentum to the Feynman parameter 
representation of the Feynman diagrams. When 
acting on Feynman diagrams one finds for exam¬ 
ple that 


d n k k-V-k _ 1 rn+2 

iW* qf.-.ql 2 lN 


(20) 


The reduction formulas necessary to express these 
integrals in terms of scalar integrals in shifted di¬ 
mensions can be found in [ [|, |[ |j • After apply¬ 
ing the reduction formula the tensor coefficient 
of each Feynman graph is expressed in terms of 
scalar integrals with trivial numerators in shifted 
dimensions similar to the approach of [ 13- 

Bases of scalar integrals 

After having reduced the scalar integrals with 
non-trivial numerators each graph is written as 
a linear combination of scalar 2- to 6-point in¬ 
tegrals in dimensions that are ranging from n 
to n + 4. By applying the scalar reduction for¬ 
mula Eq. © and explicit formulas for (n + 2 m)- 
dimensional triangle functions with light-like legs 
these integrals can be mapped to 

jn T7i jn rn+2 rn+2 ( r ) c X\ 

-*2 5 1 3,l$z2mass i ^3,3mass ’ -*4 > -*5 \^J 

Now each Feynman diagram is decomposed into 
a linear combination of these functions. As each 
of the graphs is IR/UV finite, all IR/UV poles 
have to vanish. It turns out that the coefficients 
of T 2 , 4 n ,i& 2 ma SS . ^Lss cancel as a whole in 
the amplitude. As we have proven in [ || that all 
higher dimensional integrals with more than four 
external legs have to vanish also the I£ +2 part has 
to vanish on a graph by graph bases. Note that 
the sorting into different scalar integrals is also a 
sorting into different analyticity classes, as differ¬ 
ent scalar integrals have different cut structures. 
No cancellations are to be expected between these 
classes and indeed are not observed. 

Result 

With the procedure outlined above all tensor 
coefficients were calculated. The results for these 
are lengthy and not presentable. When combin¬ 
ing them to form the actual H—f and H— helicity 
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amplitudes many interesting cancellations were 
observed which are due to the same relations that 
led to a compact result in the 6 -scalar case. 

We find for the ++ amplitude 

M ++ = £1—y 

312 4ft 

X j^i'l (S-7T47T5 ) S 7r57r6 5 S 127T3 ) 

+F\ (Sl7i- 3 ,527T 3 5 S127T3 ) 

2 (SlTJ-g J ^3 7l - 4 ,5i 7l - 37r4 ) 

+2 -F2(si27T 3 ,527T4 1T5I S17T3J S 7r47r 5 ) 

^‘2 (^1773^4 7 S27r 4 7T5 , ^1X3 : ^27T5 ) 

~F 2 (si 

7 T 37 T 4 ? ^27T37T4 ? ^77377 4 ? ^>775776 ) 

+2F 2 (si 

7T37T4 ? ^27T37T4 5 ^ 3^4 5 ^775 776 ) 

+ (1 <-+ 2) | (24) 

with 


|c(l, 7T3, 2, 7T4, 7T5, 7T6) 


^(1)2,3,4,5,6) = 2 

+ tr+(3456)R 5 + tr+(1234)Bi 

S34 

tr+(6123)R 2 +tr+(4561)R 4 

Si 6 

and 


^i(si,s 2 ,s 3 ) = - log 




F 2 (s 1, s 2 , S3,S4) 




—Li 2 



SlS2\ 

S3S4J 


(25) 


(26) 


Here tr + (j ...) = tr((l+ 75 ) .. .)/2. The ampli¬ 

tude has indeed a simple representation in terms 
of reduction coefficients as conjectured above. 


4. Conclusion 

It was pointed out that NLO corrections to 
multi-particle/jet production at future colliders 
should be included to have reliable theoretical 
predictions. On the way to understand the orga¬ 
nization of 6 -point processes at the one-loop level 
analytical results have been presented for 6 -scalar 
scattering in the Yukawa model. Very compact 


expressions could be obtained by understanding 
the cancellation mechanisms of scalar integrals. 
In the case of a nontrivial tensor structure two- 
photon scattering in the gauged Yukawa model 
has been investigated. A result has been obtained 
in terms of tensor coefficients by exploiting the 
Lorentz, gauge and analyticity structure as far as 
possible. For the H—b helicity amplitude a very 
compact representation could be obtained. The 
H— amplitude is presently under investigation. 

These results are steps towards a NLO descrip¬ 
tion of 2 —> 4 processes at future colliders. For a 
realistic QCD amplitude a lot of work remains to 
be done. 
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